UlH#.nLc.  uuK 


THE  EFFECT  OF  INTERNAL  RELAXATION  ON 
OPTOACOUSTIC  CONVERSION  IN  LIQUIDS 

Charles  Thompson  and  Henry  E.  Bass 

Physical  Acoustics  Research  Laboratory 
University  of  Mississippi 
University,  Mississippi  38677 
PARGUM  88-01 


DTIC 

ELECTEI 
FEB  2  2 1983 1 


% 


THE  UNIVERSITY  OF  MISSISSIPPI 

PHYSICAL  ACOUSTICS  RESEARCH  GROUP 

DEPARTMENT  OF  PHYSICS  AND  ASTRONOMY 


>ISTHlBUTION  STATEMOlTii 
Approrad  for  public  nlsccMV 


88  2  22  0  « 


THE  EFFECT  OF  INTERNAL  RELAXATION  ON 
OPTOACOUSTIC  CONVERSION  IN  LIQUIDS 


Charles  Thompson  and  Henry  E.  Bass 

Physical  Acoustics  Research  Laboratory 
University  of  Mississippi 
University,  Mississippi  38677 
PARGUM  88-01 


5  February  1988 

Technical  Report 
ONR  Contract  N00014-84-C-01 93 


Approved  for  public  release: 
distribution  unlimited 


Prepared  for: 

OFFICE  OF  NAVAL  RESEARCH 
DEPARTMENT  OF  THE  NAVY 
ARLINGTON,  VIRGINIA  22217 


DTIC 


SELECTE 
FEB  2  2  1988 


H 


DISTRIBUnON  STATEMENt  A 

Approved  for  public  release} 
DlstTlbutlon  Unlimited 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704  0188 


la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


1b  RESTRICTIVE  MARKINGS 


3  OISTRI8UTION/AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distributior 
uni imited 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


PARGUM  88-01 


6b  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 
(If  applicable) 

Office  of  Naval  Research 


6c.  ADDRESS  {City,  State,  and  ZIP  Code) 

Department  of  Physics  and  Astronomy 
University,  MS.  38677 


7b  ADDRESS  (C/ty,  State,  ar>d  ZIP  Code) 

Physics  Division,  Code  1112 
Arlington,  VA  22217-5000 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 


Be.  ADDRESS  (C/ty,  State,  and  ZIP  Code) 


8b  OFFICE  SYMBOL  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 

N00014-84-C-0193 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 

61153N11 


PROJECT 

NO 


(WORK  UNIT 
ACCESSION  NO 


4126936 


1 1  title  (Include  Security  Classification) 

The  Effect  of  Internal  Relaxation  on  Optoacoustic  Conversion  in  Liquids 


12  PERSONAL  AUTHOR{S) 


'  COSATI  codes  18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

"fTeLD  I  GROUP  I  SUB-GROUP 


Optoacoustlcs,  Relaxation  in  Fluids 


19  ABSTfl^CT  {Continue  on  reverse  if  necessery  er>d  identify  by  block  number) 

In  the  optoacoustic  effect,  optical  energy  is  absorbed  by  a  medium  and  converted 
to  translational  energy,  creating  an  acoustic  signal.  In  most  liquids,  the  conversion 
process  takes  place  very  rapidly  (ps  or  fs)  and  the  nonradiative  decay  time  is  much 
shorter  than  the  resulting  acoustic  signal.  Theories  have  been  developed,  based  purely 
on  geometrical  considerations,  which  accurately  predict  the  optoacoustic  signal  in  such 
cases.  The  results  show  that  for  a  short  pulse  of  optical  energy  the  time  dependence 
of  the  acoustic  pressure  is  determined  by  the  size  and  shape  of  the  excitation  zone. 

This  research  was  concerned  with  the  actual  energy  transfer  process  and  how 
variations  in  the  process  which  exist  from  one  liquid  to  another  will  affect  the  opto¬ 
acoustic  signal.  The  purpose  of  this  study  was  to  observe  the  optoacoustic  signal  in 
a  situation  where  nonradiative  decay  is  slow  enough  to  compete  with  the  geometry  of  the 
excitation  zone  in  determining  the  time  dependance  of  the  acoustic  pressure,  and  to 


20  distribution /AVAILABILITY  OF  ABSTRACT 
C3uNCLASSIFIED/UNLIMITE0  □  SAME  AS  RPT 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

L.  E.  Hargrove 


□  OTIC  USERS 


21  ABSTRACT  SECURITY  CLASSIFICATION 


22b  TELEPHONE  (/nc/ude  Area  Code)  22c  OFFICE  SYMBOL.  . 

(202)  696-4221  ONR  Code  1112 


DO  Form  1473,  JUN  86 


Previous  editions  are  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


Kt 


X 


liW 

si* 


I  j|l*.M*a*M*  A*  A*ft 


ABSTRACT 

THE  EFFECT  OF  INTERNAL  RELAXATION  ON 
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Arkansas,  1985.  M.S.,  University  of  Mississippi,  1987 

Thesis  directed  by  Professor  Henry  E.  Bass 


In  the  optoacoustic  effect,  optical  energy  is  absorbed  by 
a  medium  and  converted  to  translational  energy,  creating  an 
acoustic  signal.  In  most  liquids,  the  conversion  process 
takes  place  very  rapidly  (ps  or  fs)  and  the  nonradiative 
decay  time  is  much  shorter  than  the  resulting  acoustic 
signal.  Theories  have  been  developed,  based  purely  on 
geometrical  considerations,  which  accurately  predict  the 
optoacoustic  signal  in  such  cases.  The  results  show  that  for 
a  short  pulse  of  optical  energy  the  time  dependance  of  the 
acoustic  pressure  is  determined  by  the  size  and  shape  of  the 
excitation  zone. 

This  research  was  concerned  with  the  actual  energy 

transfer  process  and  how  variations  in  the  process  which 
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I.  Introduction 

The  generation  of  acoustical  energy  by  an  optical  energy 
source  was  observed  by  Alexander  Graham  Bell  over  one  hundred 
years  ago .  ^  Since  that  time  the  optoacoustic  effect  has  been 
used  in  investigations  of  optical  spectroscopy^ ,  acoustic 
spectroscopy^,  and  molecular  energy  transfer.^  One  advantage 
of  optoacoustic  generation  is  the  lack  of  physical  presence 
of  a  transducer  in  the  medium.  This  property  can  be  combined 
with  optical  sensing  techniques  such  as  monitoring  probe  beam 
deflection  to  produce  a  noncontact  system  for  producing  and 
observing  acoustic  signals.^ 

The  explanation  of  the  optoacoustic  effect  in  liquids  on 
a  molecular  level  begins  with  the  absorption  of  photons  by 
the  molecules  of  a  fluid.  Photon  absorption  leaves  a 
molecule  in  an  electronically  excited  state.  Dissipation  of 
this  excitation  energy  can  be  accomplished  by  several 
mechanisms  including  nonradiative  decay,  radiative  decay  and 
photochemical  processes.  In  nonradiative  decay,  interactions 
between  molecules  transfer  the  excitation  energy  to 
rotational,  vibrational,  or  translational  modes  of  the 
molecules  at  a  finite  rate.  Energy  in  rotational  or 
vibrational  modes  can  in  turn  decay  to  translation.  The  net 
increase  in  translational  energy  produces  a  temperature  rise 
in  the  fluid.  If  the  radiation  source  is  pulsed,  the  fluid 


will  cool  following  each  pulse.  The  pressure  variation  which 
accompanies  these  heating  and  cooling  cycles  is  observed  as  a 
series  of  acoustic  pulses. 

A  finite  decay  time  arises  from  relaxation  effects  in 
the  energy  transfer  process.  At  equilibrium,  the  total 
energy  of  the  fluid  molecules  is  distributed  among  their 
external  (translational)  degrees  of  freedom  and  their 
internal  degrees  of  freedom:  electronic,  vibrational  and 
rotational  modes.  When  optical  energy  is  absorbed  by 
electronic  states  some  of  the  energy  must  flow  to  the  other 
degrees  of  freedom  in  order  to  restore  equilibrium.  The 
structure  of  the  molecules  determines  the  energy  levels 
present  and  how  efficiently  molecular  interactions  transfer 
energy.  The  time  to  reestablish  equilibrium  therefore  varies 
from  one  substance  to  another. 

Each  energy  exchange  process  between  two  molecular 
degrees  of  freedom  is  described  by  a  relaxation  equation  of 
the  form: 


f  =  i(E-zJ 


where  E  is  the  value  of  the  energy  in  one  degree  of  freedom, 
Eq  is  the  value  E  would  have  at  equilibrium,  and  X  is  the 
relaxation  time  (defined  as  the  time  required  for  E  to  reach 
the  value  E^ (exp  -1))  which  characterizes  the  process.  This 


equation  also  describes  the  temperature  associated  with  each 
degree  of  freedom. 


We  can  write  the  specific  heat  of  the  fluid  as 

^p  ~  ^int  Cext  Cint  ~ 

where  is  the  contribution  to  the  specific  heat  from  a 
particular  internal  degree  of  freedom.  If  we  assume  small 
temperature  variation  we  can  write 

E  ~  E  o  —  C  ext  ~  c) 


where  T  is  the  temperature  associated  with  translation  and  Tq 
is  the  equilibrium  temperature.  For  the  external 
temperature,  then,  we  have 


dT 

dt 


=  i  (T  -  Tcj  . 
X 


Depending  on  the  time  scale  on  which  observations  are 
made,  not  all  of  the  relaxation  processes  will  be  important. 
For  example,  studies  of  excess  ultrasonic  absorption  have 
shown  that  vibrational  relaxation  times  are  generally  greater 
than  rotational  by  a  factor  of  about  100.^  This  means  that 
over  the  period  of  time  in  which  rotational  relaxation 
occurs,  little  energy  is  exchanged  by  the  vibrational  modes 
and  their  contribution  to  the  specific  heat  can  be  ignored. 

On  a  longer  time  scale,  where  vibrational  relaxation  is 
important,  rotational  modes  establish  equilibrium  quickly 


enough  to  be  considered  instantaneous  so  that  they  can  be 
grouped  with  translation  as  an  external  degree  of  freedom."^ 

Existing  treatments®' ^  of  the  optoacoustic  effect  in 
liquids  have  assumed  time  scales  long  enough  that  none  of  the 
relaxation  processes  in  nonradiative  decay  will  be  observed; 
the  transition  from  optical  to  translational  energy  is 
instantaneous.  In  this  case  the  temporal  profile  of  the 
observed  pressure  pulse  is  determined  by  factors  related  to 
the  pulsed  laser  used  as  the  energy  source.  In  the  theory  of 
Lai  and  Young®,  the  important  time  scales  are  the  laser  pulse 
length  Tp  and  the  acoustical  transit  time  across  the  beam 
radius  =  R/c  where  c  is  the  acoustic  velocity,  and  R  is  a 
measure  of  the  laser  beam  radius . 

The  purpose  of  this  study  is  to  examine  the  optoacoustic 
effect  on  a  time  scale  short  enough  to  allow  observation  of 
internal  relaxation  effects.  This  requires  that  the 
relaxation  effects  take  place  over  times  greater  than  both  l, 
and  fp.  The  two  liquids  used  are  propanol  and  carbon 
disulfide  (CS2) •  They  were  chosen  because  of  the  large 
difference  in  vibrational  relaxation  times.  Ultrasonic 
absorption  measurements  have  shown  that  CS2  has  a  vibrational 
relaxation  time  much  longer  than  most  liquids. This  makes 
it  possible  to  observe  both  optoacoustic  signals  where 
relaxation  affects  the  profile  and  those  where  it  does  not 
with  the  same  experimental  apparatus.  Instead  of  changing 


the  time  scale  of  the  experiment,  we  change  the  time  scale  of 
the  relaxation  by  using  two  liquids. 

Section  II  describes  the  theory  behind  the  experimental 
measurements  and  the  development  of  a  theoretical  expression 
for  the  optoacoustic  effect  in  CS2- 

Section  III  describes  the  experimental  apparatus  and  the 
results  of  the  measurements. 

Section  IV  is  involved  with  a  numerical  evaluation  of 
the  expression  obtained  in  Section  II. 

Section  V  compares  the  results  of  the  experiment  to  the 
numerical  calculations. 


II.  Theory 


A.  Phenomenological  approach 

Prediction  of  the  optoacoustic  signal  in  CS2  can  be 
approached  from  a  phenomenological  point  of  view  by  comparing 
the  energy  transfer  process  to  that  of  a  liquid  such  as 
propanol  where  relaxation  takes  place  very  quickly.  The 
observed  signal  in  propanol  is  very  similar  to  the  signal 
predicted  by  Lai  and  Young.®  This  indicates  that  optical 
energy  is  converted  to  translational  quickly,  so  that  the 
pressure  gradients  observed  are  due  primarily  to  the  shape  of 
the  source.  This  is  supported  by  an  estimated  value  of  the 
vibrational  relaxation  time  for  propanol  (0.02  ps)l^  which  is 
much  shorter  than  or  Xp  for  the  present  experiment .  For  a 
similarly  shaped  source,  we  can  expect  the  initial  part  of 
the  signal  to  be  the  same  in  CS2  because  the  relaxation 
process  begins  immediately  upon  absorption.  It  is  during  the 
early  part  of  the  process,  when  the  excited  states  are  highly 
populated,  that  the  greatest  amount  of  energy  is  transferred. 
At  some  point,  however,  the  signals  will  begin  to  differ  as 
translational  energy  continues  to  be  produced  in  CS2  after 
the  populations  of  excited  states  in  propanol  have 
essentially  reached  equilibrium. 

Consider  Eq.  (4)  for  the  external  temperature  T. 

Suppose  that  the  internal  temperature  T'  is  raised  suddenly 
from  Tq  to  Tj^.  The  equilibrium  value  for  T  is  also  raised; 
therefore,  energy  must  flow  into  the  external  degrees  of 


6 


freedom  to  restore  equilibrium.  The  energy  can  only  come 
from  the  internal  degrees  of  freedom,  however,  which  means 
they  do  not  remain  at  .  When  equilibrium  is  restored  it 
will  be  at  a  temperature  T2  between  and  . 

Conservation  of  energy  allows  us  to  write^^ 


Cj,nt  (l"!  To)  ~  C  int  (T  '  TJ  +  C ext  (T~To| 


At  eguili±»rium. 


^int  (T  1  T'j  —  Cext  (T~To) 


^int  (T  1  T;)  Cext  (T2~To)  ^ 


°  c  ""re  ^  c  "^rc 

int  ^  '-ext/  \'-int  ^  "-ext 


which  gives 


rr.  ™  Cint  .  _  Cext 

T2  =  Ti  +  To 


The  differential  equation  for  T  is 


=  -  (t-T’)  . 
dt  X 


Adding  Cint"^  to  both  sides  of  Eq.  (5)  we  have 


-  C;,,T'  =  CoT  -  CextT. 


(I 


Substituting  this  in  Eq.  (7)  we  have 


dt  1 


The  observed  relaxation  time,  thus,  is  not  l  but 


X.  -ii: 


With  the  use  of  Eq.  (3)  we  have 


-§  =i(E-E,) 

dt  X' 


which  has  the  solution 


E  =  E- 


dE  _  ^26 
dt  - 


This  means  that  the  translational  energy  which  appears  per 


unit  time  decays  exponentially  with  relaxation  time  x'  . 
Experimentally  then,  we  should  expect  to  see  an  exponential 
decay  superimposed  on  a  signal  due  to  the  shape  of  the 
source . 

The  above  assumes  that  all  of  the  internal  degrees  of 
freedom  have  the  same  relaxation  time.  In  real  systems, 
there  are  multiple  processes  taking  place  simultaneously. 
The  next  section  discusses  a  method  to  treat  multilevel 
systems . 

B.  Relaxation  Equation 

In  1972  Bauer^^  introduced  a  theory  accounting  for 
energy  transfer  processes  which  affect  the  optoacoustic 
signal  in  multilevel  systems.  The  matrix  notation  involved 
is  well  suited  to  numerical  calculations  and  can  be  used  to 
describe  the  effect  of  an  arbitrary  number  of  optical  and 
collisional  transitions.  His  formulation  assumes  gas  phase 
interactions  buc  is  applicable  to  liquid  with  minor  changes 
The  collisional  transitions  which  result  in  energy  transfer 
are  of  the  form 


Mi  +  Mj  M)^  Mi 


where  i...l  include  any  energy  level  of  the  fluid,  and  c 
kpj'  are  the  forward  and  backward  rates  of  reaction  Ct.  Bauer 
introduces  a  set  of  progress  variables  so  that 

n  (kaXiXj  -  ka’XkXi) 


where  n  is  the  total  number  of  molecules  and  are  the 

mole  fractions.  This  equation  has  been  modified  to  eliminate 
the  pressure  p  which  is  applicable  only  for  gases. 

Assuming  small  perturbations,  is  expanded  about 
equilibrium.  A  matrix  notation  is  then  introduced  which 
allows  summation  of  all  the  coupled  reaction  equations  into  a 
single  equation  for  the  collisional  change  in  population 
number,  ncoii.  Since  the  change  in  population  is  a  measure  of 
how  much  energy  is  exchanged  between  the  internal  and 
external  degrees  of  freedom,  ncoii  can  be  related  to  the 
temperature  variation.  The  matrix  equation  is  then  solved 
for  the  temperature  variation  as  a  sum  of  terms,  each  with  a 
particular  relaxation  time,  t,  and  amplitude,  A. 

CS2  has  been  modeled  as  a  three  level  system,  as  shown 
in  Figure  1.  The  highest  level  is  the  electronic  state  which 
is  excited  optically.  It  is  given  an  energy  of  2.96  X 
10^  cm“l  above  the  ground  state  to  correspond  to  the  frequency 
of  the  nitrogen  laser  radiation.  Between  this  level  and  the 
ground  state  is  included  a  vibrational  mode  which  is  doubly 
degenerate.  Its  energy  is  397  cm“^.^^  Although  there  are 
other  vibrational  modes  present  in  CS2r  Andreae  et .  al.^® 
have  shown  that  relaxation  of  the  vibrational  specific  heat 
of  CS2  can  be  adequately  described  by  a  single  relaxation 
time  which  they  measured  to  be  2.04  ns.  The  rate  k^o  has 
therefore  been  given  the  value  that  they  measured  and  the 
other  modes  have  been  neglected.  The  other  two  rates,  k2o 
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)c2i,  will  be  varied  later  to  give  the  best  fit  to  experimental 
data . 

With  a  three  level  system  and  one  conservation  law 
(conservation  of  total  number  of  molecules)  a  double 
relaxation  process  is  expected.  Solution  of  the  coupled 
reaction  equations  results  in  two  relaxation  times,  and  X2, 
and  two  relaxation  strengths,  and  A2 ■  Calculation  of  these 
values  was  accomplished  with  the  use  of  a  program  written  by 
Manaf  Ali  in  1985  for  his  dissertation.^^  The  program  was 
originally  compiled  to  calculate  pulsed  spectrophone  response 
in  sulfur  heixaf luoride  using  the  Bauer  matrix  notation.  It 
was  therefore  easily  modified  to  the  present  application. 

C .  Wave  Equation 

According  to  Hutchins  and  Tam^,  we  can  write  the 
inhomogeneous  wave  equation  for  the  optoacoustic  pressure  as 


1  3 


3t‘ 


-  V 


P  = 


K  +  K 

‘'“at  .  2 


at' 


/ 


(14) 


where  p  is  the  acoustic  pressure,  I  the  laser  intensity,  and 
c  the  acoustic  velocity.  The  first  term  on  the  right  hand 
side  of  Eq.  (14)  is  due  to  thermal  expansion  where 


K.  = 


ap 


(15) 


a  is  the  optical  absorption  coefficient,  P  is  the  volume 


'A 


I 


A 


*-0 


A 


13 


expansivity  and  Cp  is  the  specific  heat  at  constant  pressure. 
The  second  term  is  due  to  electrostriction  where 


Kp  = 


2ncLC 


Y  is  the  electrostrictive  coefficient,  n  is  tiie  index  of 
refraction,  and  c-  is  the  speed  of  light  in  a  vacuum. 

If  we  introduce  a  potential  function,  ())(r,t),  that 
satisfies 


^  -  V  <{)  (r,t)  =  I  ir,t) 

c  at 


then  we  can  write  the  acoustic  pressure  as  a  sum  of  thermal 
expansion  and  electrostrictive  terms  given  by 


3* 

P.  -  K,^ 


Pe  -  Ke 


^  _  ill  ^Pa 

.2  '  Ka  5t 


In  the  present  case  the  ratio  of  Kr  to  is  small  so  that  we 
will  ignore  the  electrostrictive  pressure  component. 

Using  the  laser  intensity  as  the  source  function  in  the 
wave  equation  ignores  the  effects  of  internal  relaxation  in 


& 


optoacoustic  conversion.  It  is  instructional,  however,  to 
note  results  obtained  with  this  procedure.  Lai  and  Young® 
have  investigated  cases  of  weak,  absorption  where  the  spatial 
length  of  the  source  is  much  longer  than  the  perpendicular 
distance  to  the  observation  point  so  that  the  problem  is 
essentially  two  dimensional.  The  laser  pulse  intensity  can 
be  considered  as  a  product  of  radial  and  spatial  profiles 
characterized  by  and  tp  respectively  as  defined  in  the 
introduction.  If  both  profiles  are  given  Gaussian  forms  and 
Tp  «  X^,  as  in  the  present  case,  then  the  resulting  pressure 
profile  is  relatively  insensitive  to  Xp .  The  rise  and  decay 
of  the  potential  function  and  therefore  the  peaks  in  pressure 
{d<^/dt)  are  governed  primarily  by  the  spatial  profile  of  the 
absorbing  region. 

We  have  modeled  the  source  term  for  the  present 
experiment  by  replacing  the  time  dependence  of  the  laser 
pulse  with  the  term 


A-  -l/t.  A'  '-/t; 
— ^  e  ■  +  — -  e 


where  A2,  A2  are  the  relaxation  strengths  and  Xj.Xt  are  the 
relaxation  times  discussed  in  the  previous  section.  With  a 
Gaussian  radial  profile  of  radius  a,  and  the  pulse 
propagating  along  the  z-axis,  we  now  write  the  source  term  as 


I  (r,  t|  = 


-./T.  A 

e  ■  +  ■ — 
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which  accounts  for  internal  relaxation  as  well  as  a  finite 
absorption  coefficient.  In  effect,  this  gives  the  laser 
pulse  a  flat  temporal  profile  (an  instantaneous  rise  time  at 
t=0) .  The  rise  in  the  potential  will  be  governed  by  the 
spatial  profile  as  it  would  in  the  absence  of  relaxation, 
since  tp  «  .  The  decay,  however,  will  depend  on  the 

relationship  of  to  and  X2  *  With  the  Gaussian  radial 
distribution  defined  this  way,  the  R  defined  by  Lai  and 
YoungS  is  equal  to  a/VT.  Xa  is  therefore  equal  to  a/V7c.  In 

this  form  ai  is  the  translational  energy  per  unit  time  per 
unit  volume  which  appears  in  the  liquid.  Combining  Eq.  (17) 
and  (21)  we  have 


1  9  ^^1  .  Ai  “t/Ti  A2  -r  /a  -az 

— r  — -  -  V  <|)*  —  e  +  —  e  e  e 

c  9t  /  Hi  h 


The  general  form  of  this  equation  is 


1  2 

_  — _  _  V  y  =  4Kq  (r,t) 
c  at" 


where  q(r,t)  is  the  source  strength. 

The  solution,  using  Green's  function,  is^^ 


q  (r^,  t  -  R/C 


where  I |  =  (ro2+z^2)l/2  ^3  ^he  position  of  the  source  point, 
R  =  I t-Tq I  ,  and  the  integrated  volume  includes  the  spatial 
extent  of  the  source  as  well  as  any  image  sources  necessary 


to  meet  boundary  conditions. 

In  this  experiment,  the  test  cell  containing  the  liquid 
is  large  enough  that  no  reflections  from  the  walls  are 
observed  until  well  after  the  original  signal.  We  therefore 
ignore  boundary  effects  except  for  the  surface  through  which 
the  laser  pulse  enters  the  fluid.  Boundary  conditions  at 
this  surface  are  met  by  including  a  mirror  image  of  the  real 
source  at  negative  values  of  Zq. 

Comparing  Eg.  (17)  and  (23) 

q  (r,t)  =  ^  I  (r,t)  . 

(25 

Inserting  this  and  Eq.  (21)  into  Eg.  (24),  we  have 


:,t|  =  ^  f 

1/  it1  J 


2  2 

e  -<t-  r/c)/i,  ,  Aa  -(r-  r/c)/Tj\ 

dVo  -  I —  e  —  e 


where 


<t)  (r,  t)  =  ^  g(r4  f  {t-  R/C) 


g(ro)  =  exp  |-r^a^  -  az^ 


is  the  spatial  profile  function,  and 


^  At  -t/T,  Aj  -t/T, 

f  (t)  =  —  e  *  +  —  e 
I,  T- 


is  the  temporal  profile  function- 
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III.  Experiment 

The  experimental  apparatus  used  in  this  study  consists 
basically  of  a  pulsed  nitrogen  (N2)  laser  used  as  the  energy 
source  to  create  optoacoustic  pulses,  a  continuous  wave 
Helium-Neon  laser  which  probes  the  refractive  index  gradient 
produced  by  the  acoustic  pulse,  and  a  photodetector  to 
monitor  the  probe  beam.  These  optical  components  are  mounted 
on  a  Newport  Research  Corporation  table  top  and  vibration 
isolation  system  as  shown  in  Figure  2.  The  N2  laser  is  a 
model  LNIOOO  atmospheric  pressure  laser  manufactured  by  PRA, 
Inc.  It  produces  pulses  of  1  mj  energy  and  800  ps  duration 
for  a  peak  power  of  1.25  MW  at  337  nm  wavelength.  The  pulse 
is  collimated  and  then  focused  into  a  quartz  fluorimeter 
cell,  which  contains  the  liquid  sample,  mounted  on  a 
translation  stage.  The  lenses  used  in  the  laser  path  axe 
fused  silica,  and  the  mirror  is  a  silver  coated  front  surface 
reflector.  The  laser  pulse's  original  cross  section  is 
rectangular,  3  mm  X  6  mm.  Its  size  as  it  enters  the  test 
cell  is  controlled  by  the  focal  length  of  the  focusing  lens. 
The  maximum  energy  entering  the  test  cell  per  pulse  is  about 
0.36  mJ.  The  pulse  energy  can  be  varied  by  placing  glass 
microslides  in  the  laser  path.  Each  glass  slide  absorbs  a 
small  percentage  of  the  energy  in  the  pulse.  Most  of  the 
data  taken  in  CS2  is  at  a  pulse  energy  of  about  26.3  |lJ. 
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These  values  were  measured  with  a  Molectron  J3-09  Joule  meter 
and  JDIOOO  Joule  meter  display. 

The  laser  pulse  is  absorbed  by  the  liquid  in  accordance 
with  Beer's  law  and  produces  an  outgoing  acoustic  pulse.  The 
acoustic  signal  is  observed  with  a  probe  beam  deflection 
technique  as  described  by  Sullivan  and  Tam.^"^  The  probe 
laser  is  a  Spectra  Physics  model  106-1  Ke-Ne  laser  of 
10.0  mW  output  power  at  632.8  nm.  The  beam  is  spatially 
filtered  and  then  collimated  at  a  diameter  of  approximately 
15  mm.  It  is  then  focused  into  the  test  cell  by  a  lens  of 
350  mm  focal  length  producing  a  beam,  waist  of  about  18  |lm.  A 
silver  coated  prism  reflects  the  converging  beam  into  the 
test  cell  so  that  it  propagates  above  and  parallel  to  the 
ultraviolet  pulse  as  they  enter  the  cell. 

The  pressure  gradient  in  the  acoustic  pulse  is 
accompanied  by  a  gradient  in  the  index  of  refraction.  As  the 
acoustic  pulse  crosses  the  probe  bea.m,  the  changing  index  of 
refraction  causes  a  deflection  of  the  beam  which  is 
proportional  to  the  gradient  of  the  acoustic  pressure 


r(r,t)  .  ^  a 

p.T  dr  or 


where  r(r,t)  is  the  angular  deflection,  n^  is  the  normal 


index  of  refraction,  n(r,t)  is  the  changing  index  of 


refraction  which  is  proportional  to  P{r,t),  the  acoustic 


pressure.  Over  a  length  of  the  probe  beam  d,  the  cumulative 
angular  deflection  is 


a  I  p(r,t)  dz 
a  ^  j  p(r,  t)  dz  . 


The  probe  beam  continues  through  the  test  cell  and  through  a 
lens  which  is  used  to  control  the  diameter  of  the  beam  at  the 
position  of  the  photodiode.  An  interference  filter  is  used 
to  block  any  scattered  light  from  the  pulsed  laser.  The 
photodiode  is  placed  so  that  any  beam  deflection  changes  its 
position  in  the  beam's  radial  intensity  distribution,  which 
is  assumed  to  be  Gaussian.  Positioning  the  photodiode  on  the 
steepest  part  of  the  Gaussian  distribution  assures  that  the 
intensity  change  is  approximately  linear  with  beam  deflection 
and  produces  a  maximum  change  in  intensity  for  a  given 
angular  deflection. 

The  output  from  the  photodiode,  S(t),  can  be  related  to 
the  angular  deflection  of  the  probe  beam  by 


S(r,t)  =  Glp(rd)  L 


where  G  is  a  constant  which  depends  on  the  sensitivity  of  the 
photodiode.  Ip'  (r^^)  is  the  lateral  spatial  derivative  of  the 
intensity  at  r^,  the  photodiode  position,  and  L  is  the 


distance  from  the  test  cell  to  the  photodiode.  From 
equations  (31)  and  (32)  we  have 


S  (r,t)  a  ^  J  P (^/t)  dz. 


Therefore,  when  the  signal  is  displayed  on  an  oscilloscope, 
we  observe  the  gradient  of  the  acoustic  pressure.  In  terms 
of  the  potential  function  defined  in  section  II 


S(r,t)  “  ^  ^  I  <I)(r,t)  dz. 


The  signal  is  digitized  with  a  Telctronix  7854  storage 
oscilloscope  which  is  triggered  by  a  signal  from  a  second 
photodiode  upon  which  a  small  portion  of  the  laser  pulse  is 
directed.  The  oscilloscope  is  interfaced  with  a  Digital 
Equipment  Corporation  Minc-23  mini  computer  where  the  data  is 
stored  and  analyzed. 

To  observe  internal  relaxation  effects  it  is  important 
that  the  pulse  cross-section  be  as  sm.all  as  possible  as  it 
enters  the  test  cell  so  that  is  smaller  than  and  T2  •  A 
short  focal  length  lens  is  therefore  desirable  to  focus  the 
laser  pulse  into  the  cell.  The  shortest  that  could  be  used 
with  this  arrangement  was  15  rrjri  because  of  the  prism 
necessary  to  reflect  the  probe  beam  into  the  cell. 
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Since  propanol  is  transparent  to  the  nitrogen  laser 
pulse,  it  was  necessary  to  add  a  dye  which  would  absorb  the 
ultraviolet  and  still  permit  passage  of  the  Helium-Neon  probe 
beam.  Enough  dye  was  added  to  make  the  penetration  depth 
similar  to  that  of  CS2  which  absorbs  strongly  at  337  nm. 

With  similar  penetration  depths  the  only  difference  between 
the  signals  observed  in  the  two  liquids  should  be  due  to  the 
difference  in  sound  speed  and  internal  relaxations. 

Typical  waveforms  in  both  liquids  are  shown  in  Fig.  3. 

To  account  for  the  differences  in  sound  speed  we  compare  the 
two  signals  on  a  scale  of  ct .  In  Fig.  3  ct  =  0.0  is  the  time 
at  which  the  oscilloscope  was  triggered.  As  expected,  the 
early  portion  of  the  signals  are  similar.  Following  this 
initial  rise,  however,  the  CS2  signal  decays  more  slowly  and 
reaches  its  minimum  at  a  slightly  later  time.  The  pressure 
gradient  in  CS2  xs  smaller  in  magnitude  than  in  propanol, 
most  of  the  time- 


I 


■ « ■ .  • « -  •  ■ .  • .  *.«  Vv","  V  W 
.  W.V'.V.S’.V.N  .NwV  .S -N 


IV.  Numerical  Calculations 
A.  Optoacoustic  Signal 

From  equation  (34)  the  expression  we  compare  to  the 
experimental  signal  is 


A  ^  ir.tl  dz  . 

If  we  define 


(35) 


O  (r,  t)  =  [  (Hr,  t)  dz 
and  use  equation  (27) ,  then 


(36) 


^1  \  f'*  1  f  <1  (ro)  f  (t  -  R/c) 

0(r,t)  =  J  dz  cVo  i ^  - 

•'0  /  \2  j  R 

(471) 


If  we  define  a  function  u(r,t)  by 


ur 


,  t)  «  -J—  f  <s/„  2!^  5(t  -  R/cl 


and  define  the  Dirac  response  of  the  probe  beamJOy 


(37) 


(36) 


U 


(r,  t)  =  f  u(r,  t)  dz 

J  n 


then 


(39) 


t-t’)  f(t')  dt'  . 


(40) 
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Writing  O  in  terms  of  U  saves  time  in  the  numerical 
evaluation  of  <1>  because  the  spatial  integration  need  only  be 
done  once  rather  than  at  each  point  in  time.  u(r,t) 
represents  the  potential  due  to  a  source  with  spatial 
dependance  g (r^)  ,  which  exists  only  at  t=0 .  The  integral  of 
u(r,t)  over  z  represents  the  effect  upon  a  length  d  of  the 
probe  beam.  The  convolution  of  the  resulting  function  with 
the  actual  time  dependence  of  the  source  produces  the 
potential  at  any  later  point  in  time. 


Evaluation  of  u(r,t)  is  accomplished  by  integrating  over 
all  space  in  a  spherical  coordinate  system  with  its  origin  at 
the  observation  point.  Defining  dVo  by 

dVo  =  R^sin0  dR  d0  d<>  =  R^dR  dQ 


we  have 


u(r,  t)  =  — - —  I  j  R  dR  dQg(ro)  6(t  -  R/c)  . 


It  is  possible  to  show  that 


5(t  -  R/c)  =  c6(ct  -  r| 


therefore , 


u(r,  t|  = 


)47t) 


R  dR  dng(rj  6|ct  -  r) 


The  integrand  is  zero  except  when  R  =  ct,  therefore 


u(r,  t)  =  -  I  i  g(rc)  sin  6  d0  ci(}) 

(4:1)'  “  “ 


Defining  X  =  cos0,  dx  =  -sin0d0, 

u(r,t|  =  -5-5-  I  I  gjrj  dx  dO . 

A  point  in  time  t  defines  a  sphere  of  radius  R=ct  about  the 
observation  point  irj  =  +  u{r,t)  is  the  sum  of  the 

contributions  from  all  source  points  which  lie  on  the  surface 
of  this  sphere.  In  the  numerical  calculations,  knowledge  of 
the  source  has  been  used  to  m.odify  the  limits  of  integration. 
By  defining  g (r^)  to  be  zero  for  r^  >  4a  and  Zq  >  20/a, 
definite  limits  for  X  <{)  are  calculated  for  each  value  of 

r  and  t.  The  length  of  the  probe  beam  considered  was 
determined  by  calculating  the  farthest  point  along  its  length 
that  would  be  affected  by  the  source  at  the  final  time  of 
interest.  The  programs  which  calculate  U(r,t)  and  are 

located  in  the  appendices.  The  pressure  ^  ^(r.t)|  and  its 

gradient  ^  ^(r,t)  also  calculated  num.erically  by 

program, s  in  the  appendices. 
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B.  Curve  Fitting  Procedure 

The  procedure  for  producing  the  theoretical  curves  was 
to  begin  with  very  short  relaxation  times  so  that  the  signal 
is  dominated  by  the  spatial  characteristics  of  the  excitation 
beam.  These  characteristics  were  varied  until  the  result 
matched  the  observed  signal  in  propanol  which  is  known  to 
have  a  very  short  relaxation  time.  Relaxation  times  and 
strengths  were  then  calculated  using  the  known  rate  for  the 
vibration  to  translation  transition,  kj^O/  assumed  rates  of 

transition  from  the  electronic  state.  Varying  the  values  of 
k2o  c^rid  k2i  changed  the  relaxation  times  and  strengths.  These 
two  rates  were  varied  until  the  resulting  relaxation 
strengths  and  times  gave  an  optoacoustic  signal  which  most 
closely  resembled  the  observed  signal  in  CS2 • 


V.  Discussion  and  Conclusion 
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The  results  of  this  study  show  that  internal  relaxation 
can  indeed  be  important  in  determining  the  time  dependence  of 
the  optoacoustic  signal  in  liquids. 

A  good  representation  of  the  signal  observed  in  propanol 
was  obtained  by  assigning  a  rate  of  490.2  X  10^  s~^  to  each  of 
the  three  transitions.  The  resulting  relaxation  times  were 
1.02  ns  and  1.47  ns.  The  pressure  gradient  obtained  with 
these  relaxation  times  is  shown  in  Fig.  4  along  with  the 
experimental  curve.  Here  the  time  scale  is  for  the 
theoretical  curve.  The  experiemntal  curve  is  overlaid  so 
that  the  peaks  occur  at  the  same  point  in  time.  Performing 
the  calculation  with  even  shorter  relaxation  times  results  in 
the  same  theoretical  curve,  indicating  that  with  these 
relaxation  times  the  signal  is  determined  entirely  by  the 
spatial  characteristics  of  the  excitation  pulse.  The  value 
of  "a”  in  this  calculation  was  50  giving  the  radius  a 
value  of  approximately  71  |r,m  at  the  exp  (-2)  value  of  the 
intensity.  The  theory  of  Heritier^  predicts  that  the  time 
separation  of  the  positive  and  negative  pressure  peaks,  which 
is  approximately  the  time  separation  of  the  points  where  the 
gradient  goes  through  zero,  will  be^”^ 


At  =  1.66 


CO 


t  + 

p 


2c  / 
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where  (Oq  is  the  radius  at  which  the  intensity  is  down  by 
exp  (-2)  .  For  the  present  case  with  Xp  =  0.8  ns,  CDq  =  71)im  and 
c  =  1223  m/s  for  propanol.  At  should  be  68.3  ns.  The  actual 
time  separation  in  the  experimental  waveform  is  72.3  ns.  The 
discrepancy  in  the  initial  rise  of  the  curves  in  Fig.  4, 
which  begins  earlier  in  the  experimental  curve  and  shows  a 
more  gradual  change,  is  probably  due  to  the  fact  that  the 
laser  pulse  is  not  Gaussian  as  modeled.  The  roughness  of  the 
computed  waveform  can  be  traced  to  edge  effects  in  the  source 
function  caused  by  truncating  the  Gaussian  distribution  at 
Zq  =  4a  and  =  20/a.  These  effects  are  magnified  in  the 
process  of  taking  two  numerical  derivatives  of  the  calculated 
function,  O.  Both  of  these  problems  can  probably  be 
corrected  by  modifying  the  spatial  profile  function,  £q. 

(28) .  However,  the  overall  time  dependence  of  the  pulse  is 
well  represented  by  the  theory.  It  should  be  noted  that  to 
avoid  noise  during  the  experimental  signal  digitization,,  the 
band  width  of  the  oscilloscope  was  reduced  to  20  MHz.  Had  a 
similar  filter  been  used  in  the  numerical  calculations,  the 
agreement  might  be  different. 

Table  1  shows  combinations  of  values  assigned  to  the 
rates  kn^,  kj^g^  k2o  s^id  the  resulting  amplitudes  and 

relaxation  times.  Theoretical  wave  forms  calculated  using 
these  parameters  are  shown  in  Figs.  5  through  8. 
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Table  1 


Figure  Rates  A2  X2 

in  10^  s“l  ns  ns 


^21 

klO 

^20 

5 

490.2 

490.2 

490.2 

-1.26x10-5 

2.03x10-'* 

1.47 

1.02 

6 

26.86 

490.2 

16.86 

-1.26x10-5 

2.03x10-4 

1.47 

22.8 

7 

26.86 

490.2 

.5686 

-1.26x10-5 

2.03x10-4 

1.47 

36.5 

8 

.2686 

490.2 

10 . 67 

-1.26x10-5 

2.03x10-4 

1.47 

91.4 

With  the  value  of  k^Q  held  constant  to  correspond  to  the 
accepted  value  of  the  vibrational  relaxation  time  of  CS2r  the 
amplitudes  and  also  remain  constant.  Table  1  is  arranged 
to  show  rate  combinations  which  give  increasing  values  of  X2. 
These  results  are  not  unique,  however.  For  instance,  the 
last  value  of  X2,  91.4  ns,  also  results  from  setting  k2i  = 

20  = 


10.37  X  10^  s~^  and  k 


.5686  X  106  s-1 


Figure  9  compares  the  pressure  gradients  from  Figures  6, 
7,  and  8  and  clearly  shows  the  effect  of  increasing  X2.  With 
the  first  positive  peak  in  each  curve  normalized  to  the  same 
magnitude,  the  negative  peak  and  the  second  positive  peak 
decrease  in  amplitude  as  t2  increases  while  all  three  peaks 
broaden  in  time. 

The  best  fit  of  the  experimental  CS2  waveform  was 
obtained  with  the  rates  which  produced  the  curves  in  Fig-  8. 

A  comparison  of  the  theoretical  and  experimental  curves  is 
shown  in  Fig.  10.  As  in  Fig.  4,  the  time  scale  is  ior  the 
theoretical  curve.  The  experimental  curve  is  positioned  so 
that  the  peaks  occur  at  the  same  point  in  time.  The  value  of 
T2  for  this  curve  is  91.4  ns.  Increasing  the  value  of  T2 
beyond  this  point  continues  to  decrease  the  relative 
magnitude  of  the  negative  peak,  which  the  comparison  shows  to 
be  larger  than  experimentally  observed.  However,  as  ^ 
becomes  larger,  the  initial  pressure  gradient  becomes 
negative.  From  equation  29  we  can  see  that,  although  A2  is 
positive  and  larger  in  magnitude  than  Aj,  when  »  Tj,  the 
first  term  becomes  dominant  and  f(t)  is  negarive  for  small  t. 
Therefore,  with  the  three  level  system  considered  here,  we 
are  limited  in  how  large  the  difference  in  and  T2  can  be 
and  still  produce  the  same  general  shape  of  the  CS2  signal. 


W  T<.s 


A  better  prediction  of  the  optoacoustic  signal  should  be 


obtainable  by  considering  a  more  comprehensive  reaction 
scheme  than  the  three  level  system  presented  here. 

The  agreement  is  better  for  the  early  part  of  the  pulse. 
The  same  discrepancy  in  the  initial  rise  that  appeared  in  the 
propanol  signal  occurs  here,  and  again  could  be  corrected  by 
modifying  the  spatial  profile  function. 

From  these  results  we  conclude  that  the  optoacoustic 
effect  can  be  used  to  investigate  relaxation  effects  in 
liquids.  Obviously,  to  observe  relaxation  in  most  liquids, 
which  have  much  shorter  relaxation  times  than  CS2,  it  will  be 
necessary  to  produce  much  shorter  acoustic  pulses  than  the 
apparatus  used  here  produces.  However,  the  theoretical 
calculations  show  that  the  effects  of  relaxation  become 
apparent  even  for  relaxation  times  less  than  for  the 
experiment.  Figure  6  shows  a  distinct  difference  from  Fig.  5 
in  the  latter  part  of  the  signal.  The  longest  relaxation 
time  (lowest  transition  rate)  involved  in  the  Fig.  6 
calculation  was  1/^20  ~  59.3  ns;  and  X2  for  this  curve  was 
22.8  ns,  which  is  less  than  (~31  ns) .  An  experimental 
arrangement  that  allows  better  focusing  of  the  excitation 


laser  should  therefore  increase  the  number  of  liquids  in 
which  it  is  possible  to  observe  internal  relaxation. 


,1  j 
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APPENDIX  A 

ROUTINE  FOR  CALCULATION  OF 
riRAC  RESPONSE  OF  THE  PROBE  BEAM 
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irac'  calculates  the  dirac  response  o-f  the 
obe  beam  as  an  integral  over  z. 

source  is  a  delta  -function  in  time,  o-lindricall 
shaoed  with  a  oaussian  radial  Dro-file,  exponent  i 


X  fij 


r^t’i 


ml=4 
nl=20 

r=2 . 5e “3 

*  print*,'  input  r  value  in  meters'' 

*  read*,r 

*  this  section  determines  the  -first  point  in  time  that 

*  there  is  a  contribution  to  U  <  the  shortest  distance  to 

*  a  source  point  )  and  how  big  the  time  increment  must  b« 

*  to  include  all  o-f  the  response. 


to=  <  r-ml *a)/c 

for  now  let  total =1 500ns  to  corespond  to  D-scope 
pictures  <50ns/div) 
total=  1500. 0e-9 
dt=  total/300 
t  <  0 )  =  t  o-d  t 

t  <  0  >  is  the  first  point  in  the  f  i  1  ie  '  t ' 
it  is  placed  1  dt  before  the  first  contributing 
time,  consequently  U(0)  is  always  8.0 

wr i t  e  < 1 1  , * ) U< 1  , 0 ) 
wr  i t  e  (  1  2  ,  *  >  U  <  2 , 0  ) 
write  < 1 3 , *  >  t  <  0 ) 

here  we  calculate  how  far  along  the  probe  beami  we 
must  integrate  untill  the  propagation  time  is  greater 
than  tf=to+500ns. 

z=  (  n  1 /a  1  p h a )  c  »sqr  t  (  <  500 . 0e-9 )  **2+2. 0*tcr*500 . 0e -9 > 

print*,'  zrriax= '  ,  z 

there  are  z*alpha  penetration  depths  in  this  distance 
and  we""!!  have  5  points  per  penetration  jiev  th 

p ds  =  n i n  t ( z  *a 1 p h  a ) 
dz  =  0 . 2/ a  1 ph  a 

print*,'  ',pds,'  penetration  depths  '  ,5*pds  ,  z-points 


as  we  step  through  time,  each  poin 


associated  with  the 


*  surface  of  a  sphere  of  radius  rho=ct(i). 

*  we  haue  contributions  to  phi  from  source  points  which  lie 

*  on  this  surface.  <whereuer  the  sphere  intersects  the 

*  source  distribution.) 

*  at  each  point  in  time  we  will  determine  the  range  of 

*  fee  and  xi  which  are  in  the  source  and  sum  the  contribution! 

*  from  source  pts.  in  this  range. 


do  1100  point=l,2 

r=r  +  ( po i n  t-1 ) *c*dt 
do  1010  s=0 , 5*pds 
z=s*dz 

p  r  i  n  t  *  ,  '  i  n  p  u  t  z  ' 
read* , z 

do  1000  i=l ,300 
t< i )=t(0)+i*dt 
r  h  o=  c  *  t  (  i  ) 

giuen  rho,  we  can  determine  the  limits  of  xi 


case  1 


c  ase  2 


if  ( r ho . 1 t . ( r-a*ml > )  go  to  1000 

if  < r ho . 1 t . < r + a*ml ) )  then 
X i max=l . 0 

X  i  m  i  n=  <  r -a«m 1 ) /r  h  o 
e  1  se 

X i max=( r  +  a«ml )/r  ho 
X i m i n= <  r -a*m 1 > /r  h  o 
en d  i  f 


dx i  =( X i max-x i m i n )/20 

in  order  to  do  the  calculation  at  the  center 
of  each  dx i ,  moue  the  first  xi  to  xi+  dxi/2-0 
ximin=  ximin  +  dxi/2.0 

do  900  .j=0,19 

X i =  ximin  +  J  *dx i 

given  rho  and  xi  we  can  determine  the  1  imits  on  "fe^e 

maximum  value  of  fee  is  p i /2 

this  section  determines  the  minimum  value 

t e s t 1 = < a*m 1 ) **2  -  (rho*xi)**2  -  r**2  + 

2*r  *r  ho*x i 


SC 
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test2=rho**2*( 1-x i **2) 

\i  ( test  1 . ge . test2)  then 
m  i  n-f  ee=0 . 0 
e  1  se 

m  i  n  -f  e  e=acos<  sqrt<  testl/te£t2)  > 
end  i  -f 

limits  on  -fee  are  min -fee  to  p  i /2  and  -min -fee  to 
-pi/2  contributions  from  negative  angles  will  be 
paired  with  positives  when  calculating  U. 

d-fee  =  <  p  i /2  -  m  i  n-f  ee  )/20 . 0 
rriin-fee  =  min-fee  d-fee/2 

do  300  k=0,19 

■fee  =  rriin-fee  +  k^d-fee 


now  we  have  r ho, xi, -fee  values 
calculate  ro,20,  and  9(ro,zo) 
ro  is  independant  o-f  the  sign  of  fee, 
zo  has  different  values  for  fc-e  and  -fee 

ro  =  sqrt(  <rho*xi  -  r)**2  (  r  h  o*sq  r  t  (  1 -x  i  **-2  > 

*cos ( f ee ) ) **2) 

zo  =  z  r  h o*sq r  t  <  1 -X  i  )  *s  i  n  <  f  e e ) 

if  zo  is  out  of  range  don  ■' t  allow  contribution  to  U 

go  instead  to  check  zo(-fee) 

if  < abs( zo> . gt . < n 1 /al pha) )  go  to  100 

g  =  e**''. -<  r  o/a)  -  al  pha*(  abs(  zo)  >  ) 

U  (  p  o  i  n  t  ,  i  )  =  U  p  o  i  n  t  ,  i  )  t  (  i  )  *(  c/  (4^pi)>*'  +  2  * 
g*dx i *df ee*dz 

zo  =  z  -  rho^sqr t ( 1 -X i **2) *s i n ( f ee > 

if  zo  is  out  of  range  don't  allow  contribution  to  U 

go  instead  to  get  new  fee 

if  ( abs( zo) . gt . < n 1 /al pha) )  go  to  300 

g=  e>-*'; -<  ro/a)  -  al  pha*<  abs  <  zo)  )  > 

U  <  p  o  i  n  t  ,  i  )  =  U  <  p  o  i  n  t  ,  i  )  t  (  i  )  *  <  c  /  <  4  *  p  i  )  )  )^  -^  2 
*  Q>dx i *df e e #dz 
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con  t i nue 


{800 

i  s 

the 

end 

ot 

a  -fee  1  oop 

con  t  i 
{900 

t  nue 
i  s 

the 

end 

ot 

a  X  i  1 oop  > 

con  t i nue 


{1000  completes  a  pair  o-f  U(i),t<i)  -for  one  value  o-f  z> 
continue 

{1010  c  omp  letes.  LKi),t(i)> 
continue 

{1100  ends  calculation  ot  U  -for  one  r  velue? 

{write  arrays  to  tiles} 
do  1002  i =1 , 300 

write  ( 1 1 ,*)U< 1 , i ) 
write  <  1  2  ,  * ) U ( 2 ,  i  ) 
write  (13,*)t(i> 
c on  t i nue 


APPENDIX  B 

ROUTINE  FOR  CALCULATION  OF 
POTENTIAL  FUNCTION 


'duhamel'  calculates  the  potential  using  duhamel 
principle  and  the  dirac  response  -from  'dirac' 


'dirac'  computes  a  301  point  array  for  U(i),t(i) 
ph i < i )  is  the  integral  over  tp  from  0  to  t(i)  of 
U<  t<  i  )-tp)f <  tp)dtp 


however  U  is  zero  for  t  less  than  t(0) 

or  t<i)-tp  less  than  t(0) 
or  tp  greater  than  t(i>-t(0) 
so  we  set  an  upper  limit  for  the  integral  of 
tp/dtp  =  < t( i )-t<0) )/dtp  =  i 


p  h i  ( i  )  -  potential 

U( i )  -  dirac  response 

f  -  time  dependant  profile  function 

taul,tau2  -  relaxation  times 

A1 ,A2  -  relaxation  strengths 

tp  -  integration  variable 


real  phi  <2,0:300) , U  <  2 , 0 : 30  0 ) , t  <  0 : 30  O  > 
real  dt,dtp,tp,f 
real  t au 1  , t au 2 , A 1  , A2 
integer  i,l,ul,q, point 

open  <  u n i t= I  4 , f i  1 e=  ^  p h i  1  '  ) 
open  (  un i  t=l 5 , f i 1 e=' ph i 2' ) 

data  ph i /602*0 . 0/ 

open  (uni  t= 1 1  , f i 1 e= ' U 1  '  , s  t  a  t  u  s= ' o 1 d ' > 
open  (uni  t  =  1 2 , f i  1 e= ' U2 ' , s  t  a  t  u  s= ' o 1 d ' ) 
open  (uni t=13,fi le='t^ , statu s= 'old') 
read(  1  3  ,  ♦^)  t  (  0  ) 

print*,'  input  values  for  taul,tau2' 
r  e  ad* , t  au 1  , t  au  2 

print*,'  input  values  for  A1 ,A2' 
read*,Al ,A2 
Al=-1 . 26e-5 
A2=  2.031e-4 


do  1  1=0,300 


reacl(  1  1  ,*)U<  1  ,  i  ) 
read(12,*>U<2, i ) 
continue 

dt=5.e-9 

note  that  setting  dt=5ns  is  specialized  to  form  of 
^dirac'  where  dt=5ns 

do  100  po i n  t=l  , 2 
do  50  i =0 , 300 

t<  i  )  =  t(0)+ i *dt 
ul  =  i 

do  30  q=l  ,2 

if  <q.eq.l)  then 
tau=tau 1 
A=A1 
e  1  se 

tau=tau2 
A=A2 
e  n  d  i  f 

two  loops  are  included,  one  for  fast  relaxations 
one  for  slow  relaxations 

if  ( t au . 1 e . 2«dt )  go  to  20 

continue 

dtp=dt 
tp=0 .0 

the  integration  is  done  by  the  trapezoid  rule 
f=<A/tau) 

ph i ( po i n  t , i )=ph i <  po i n  t , i ) po i n  t , i ) *f 
*dtp/2 . 0 

this  is  the  first  term  when  using  the  trapezoid  rule 
last  term  {tp=ul*dtp>  will  be  zero  because  U<step,0) 


IS  zero 


do  15  1 = 1 , u 1 
tp=l *d 

f  =  <  A/ 1  au ) *e  X  p  < - 1  p/ 1  au ) 
if  ( t p/ t au . g t . 20 . 0 )  go  to  1? 
because  f  will  be  small  for  remaining  I's 

ph  i  (  s  t  ep  ,  i  )=ph  i  <  s  t  ep  ,  i  ^  +  U s  t  ep  ,  u  1  - 1  )  *f  ^^d  tp 


con  t i nue 
go  to  30 
continue 


start  here  -for  fast  relaxations,  linear 
interpolation  is  done  between  U< i )  and  U<i+1) 
dtp=dt/25 

do  30  1  =  l  ,  u 1 

slope=<U(step,ul-l+I)-U(step,ul-l) )/25 
tq=< 1 -1 )*dt 

tp=tq 

f — < A/ tau)*exp(-tp/tau) 

phi (step, i )=phi (step, i >+U(step,ul-1 +1) *f *dtp/2 . 0 

tp=tq+25*dtp 

f=(A/tau )*exp(-tp/tau> 

phi (step, i )=phi (step, i >+U(step ,ul -1 >*f*dtp/2.0 


do  25  n=l ,24 
tp=  t  q  +  n  *d  t  p 

f=( A/tau ) *exp ( -tp/tau ) 

if  ( t p/ t au . 9 t . 20 )  go  to  38 


phi (step, i )=phi (step,! )+(LKstep,ul-l+l)-slope*n) 
*f *dtp 

con  t i nue 

continue 
continue 
cent i nue 
do  101  i =0,300 

wr  i  t  e  (  1  4  ,  * )  p  h  i  (  1  ,  i  ) 
wr  i  t  e  (  1  5  ,  *  >  p  h  i  (  2  ,  i  > 
con  t i nue 


numerically  dif+erentiates  the  potential 
to  get  the  pressure 


real  phil< 0:100), phi2(0:100>,t(0;100) 
real  presl(0:99) ,pres2(0:99),tp<0:99) 
real  delta 
integer  n 


open  (uni  t=13,ti  le='t'  , statu s=  •'old') 
open  (uni  t=14,-fi  le=''phi  I  '  , statu  s=  •'old') 
open  (un'  t=15,-fi  le='phi2'  ,status='old') 
do  900  n=0 ,100 
r e ad( 1 4 , *) ph i  1 ( n  ) 
read( 15,*)ph i 2(n) 
read( 1 3 , *) t ( n  ) 
con  t i nue 


open  (uni  t=16,'fi  le='tp'  , statu  s='  new'  > 
open  (uni  t=17,-fi  le  =  'presl'  , statu  s='  new'  ) 
open  (uni  t= 1 8 , t i 1 e= ' p r e s2 ' , s t a t u s= ' new ' > 

do  1000  n=l , 1 00 

del ta=  ( t(n)-t(n-l ) ) 

presl(n-l)=  (phi l(n)-phi l(n-l)) /de 1 t  a 
pres2(n-l )=  (phi2(n)-phi2(n-l)) /de 1 t  a 
wr i t  e ( 1 7 , * ) p  r  e  s  1  ( n - 1  ) 
wr  i  t  e  (  1  8  ,  * )  p  r  e  s 2  (  n  - 1  ) 

tp(n-l)=  t(n-l)  +  delta/2.0 
wr i t  e ( 1 6 , * ) t  p ( n - 1 ) 
con  t i nue 
end 


APPENDIX  D 


ROUTINE  FOR  CALCULATION  OF 
GRADIENT  OF  THE  ACOUSTIC  PRESSURE 


u 


■.V 


% 

Ss 


numerically  di-fferentiates  the  pressure 
to  get  it's  gradient 


real  pr e  s 1 ( 6 : 99) , pr e  s2( 0 : 99) , t p ( 0 : 99 ) 
real  dpres( 0 : 99) , de 1 ta 
i n  teqer  n 


del ta=< 1 140.0*5.e-9) 

open  (uni  t=  1  6  ,  -f  i  1  e= '  t p  '  ,  s t a t u s=  '  o  1  d '  ) 
open  (uni t= 17, til e= ' presl ' , statu s='ol d' ) 
open  (uni t=13,ti le='pres2' ,status='old') 

do  900  n=0,99 
r  e  a  d  (  1  6  ,  ^^ )  t  p  (  n  > 
r  e  a  d  (  1  7  ,  * )  p  r  e  s  1  (  n  ) 
r  e  a  d  (  1  S  ,  *  )  p  r  e  =  2  (  n  ) 
c  o  n  t  i  n  u  e 


open  (uni  t  =  1  9  ,  -f  i  1  e=  '  dp  r  es '  ,  s  t  a  t  u  s= '  new '  ) 
do  1000  n=0 , 99 

dpres<  n )  =  ( presl ( n ) -pres2( n ) )/de 1  ta 
wr i t  e ( 1 9 , * ) dp  r  e  s ( n ) 

0  continue 
end 


Tim 


